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Abstract 

In this letter we study the process of gluon fusion into a pair of Higgs 
bosons in a model with one universal extra dimension. We find that the 
contributions from the extra top quark Kaluza-Klein excitations lead to a 
Higgs pair production cross section at the LHC that can be significantly 
altered compared to the Standard Model value for small values of the 
compactification scale. 



1 Introduction 

In spite of the great experimental successes of the Standard Model (SM), we 
still do not have a direct test of its symmetry breaking sector. In fact, one may 
say that there is a tension arising from the indirect bounds on the Higgs mass 
coming from the loop structure and precision measurements, which favours a 
Higgs mass that is already excluded by direct searches. In addition, we already 
know that the SM is incomplete, since it does not provide for a viable dark 
matter candidate and for neutrino masses and mixings. 

On a more theoretical side, the SM is not satisfactory due to the triviality 
and hierarchy problems in the Higgs sector. These problems suggest that the 
Higgs sector should be viewed as an effective theory valid up to a certain energy 
scale. The traditional solutions to these problems used to be represented by two 
broad classes of models: supersymmetry and technicolor. 

Recently a third class of solutions has been proposed, involving the existence 
of extra space- like compact dimensions [T]. This class can be further divided 
in three different classes depending on the fields that can propagate in the 
extra dimensions and the geometry of these extra dimensions. In this letter we 
will concentrate on the so-called models of Universal Extra Dimensions (UED), 
where all fields can propagate in the flat compact extra dimensions f?. An 
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important property of UED arises from momentum conservation in the extra 
dimensions. This implies that KK number is conserved in all tree level vertices. 
As a consequence contributions to electroweak observables arise only from loops 
of KK particles allowing compactification scales as low as 500 GeV [3] . 

We are interested in the consequences of this model for Higgs boson produc- 
tion at hadron colliders. The analysis of single Higgs production in gluon fusion 
process was done by Petriello A] and a significant enhancement compared to 
the SM was found. 

Higgs pair production via gluon fusion in the SM was studied in ^ and 
is an interesting process since it could give information on the Higgs boson 
cubic coupling. Hence it is important to examine possible deviations from SM 
predictions in different models. For instance, Higgs pair production in Little 
Higgs models was studied in [6]. 

In this paper we study the modifications of the Higgs pair production cross 
section via gluon fusion in UED. 

2 Model and relevant masses and couplings 

In models of UED all fields are allowed to propagate in the bulk and hence they 
all have an associated Kaluza-Klein (KK) tower. We will work in the case of one 
additional compact dimension. In order to retain the zero modes corresponding 
to SM particles it is usual to compactify the extra dimension in an orbifold 
5*1/^2, defined by the identification y ~^ y + 7ri?, where y is the 5th dimension 
coordinate and R is the compactification radius, and demand that the fields 
with zero modes to be even under the transformation y — )■ —y. 

After compactification the relevant fields for our purposes (Higgs doublet H, 
top quark singlet t, top quark doublet Q and gluon field G) will have the usual 
KK expansion: 

00 

/^(a:^2/) =i?°(a;)x(°) +^[i/(")(a;)x(")(2/)] (1) 

00 

t{x^,y) = 4(x-)x(°) +Y.'^tf{x)x^-\y) + 4")(:r)0(")(y)] (2) 

00 

G^(x^^/) = G;i(")(a;)x(°) + ^ t^M^"^ (3) 

00 

Q{x^,y) ^ g°i(a;)x(°) +$^[Oi"^(x)x("Hy) +Qi"^(x)0(")(y)] (4) 

00 

GUx-,y) = J2[Gf"\x)^^-\y)] (5) 

n=l 
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where X^"'\y) and ^'"'(?/) are orthogonal basis: 

x(")(y) = ^cosf , X(°) = ^, <^("H?/) = ^Binf . (6) 

The coupHngs of KK top quarks with gluons and Higgs field that enter in 
our computation are derived from the lagrangian: 

/TT-R p 
dy / df^xiiQIpQ + itipt + [\\Qia2H*t + h.c]}, (7) 

where the covariant derivative is: 

Ip^T^^dM-ig^T'^Gl,), (8) 

M is a Lorentz index with values M = 0, 1, 2, 3 and 5 (we will use lower case 
greek index like /i = 0,1,2,3 to denote the usual non-compact dimensions), 
= (7'',i7^) are the 5-d Dirac matrices, gs and X\ are the 5-d QCD and 
Yukawa coupling constants respectively and T°' is the usual color group gener- 
ator. 

Considering only the coupling to the Higgs zero mode, which after sponta- 
neous symmetry breaking is written in terms of its vacuum expectation value 
V = 246 GeV as 

H^°' = ^,( lu..) (9) 



V2 V v + h'^'^Hx) 

A' 

and using the relation between the 5-d and 4-d top Yukawa coupling, A = 

with the zero mode top quark mass given as usual by nit = one finds that 
the mass eigenstates of the KK top quark tower are 

These mass eigenstates are related to the original states by 

Tn -Un \^^in) j -U, y ^^^^ j (11) 

where and denote the upper components of the doublets and 
respectively. The orthogonal matrices and are given by 

\ — sin ^^2" cos ^^2~ / \ sm ^^2" ~ cos ^^2~ I 

where sina^") = and cosa^"^ = ^^ii- ^ith m?„ =rni -\-mi and m„ = 

The two top KK towers mass eigenstates, f^"^ and , have a degenerate mass 
given by mt,„. 



3 



The couplings of the top KK tower states with gluons are simply given by 



oo 



(13) 



whereas the couphng to the zero mode Higgs boson can mix 4"' and 4"^: 



n=l 



+ h.c. 



(14) 

Notice that the top KK Yukawa coupHngs are proportional to the top quark 
mass and hence their effects decouple for higher KK modes. 



3 Model implementation and results 

We implemented the new particles and couplings in FeynArts [7^ for an arbi- 
trary number of KK modes. We then use FormCalc [8J to perform the com- 
putation of traces and the reduction of the tensor one-loop integrals to scalar 
Passarino-Veltman integrals [9]. Finally, LoopTools [10] computes numerically 
the integrals and CUBA [11] integrates over phase space to find the cross sec- 
tion. We verified that in the case of single Higgs production in UED, where 
only a triangle diagram contributes, the program reproduces both analytically 
and numerically the results obtained by Petriello [4] . We have also checked our 
code with the SM Higgs pair production [5] . 

In Figure [Tj we show the diagrams that are computed for one top quark 
KK level. Notice the presence of 2 top KK excitations for each level and their 
mixture through the Yukawa coupling. 

In this work we will consider compactification scales 1/i? as low as 500 GeV, 
as allowed by electroweak precision measurements [3] (see also [12] for bounds 
coming from 6 — > S7 processes) and first compute the deviations from the SM 
for a fixed partonic center-of-mass energy as a function of the Higgs boson 
mass Mh for diff'erent values of the compactification scale. We will include in 
the calculation a number n of KK levels such that to„ < 10 TeV for a given 
compactification scale, where one expects that the 4-dimensional effective theory 
starts to loose its validity [2]. For instance, we considered the contribution of 
20 KK levels for 1/i? = 500 GeV. In practice, the convergence for large values 
of n is very rapid. 

In Figure [2] we show the differences between the SM and UED contributions 
for the triangle and box diagrams separetely. We fix \/I =1.0 TeV and 1/i? = 
500 GeV for illustration. 

The triangle contribution can be understood analytically; it is given by the 
difference in the triangle amplitude: 

'^UEpjgg ^ HH) - a^^,{gg ^ HH) ^ \Asm + Akk\^ - \Asm\' 
a^,,{gg^HH) \Asm\^ 
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Figure 1: Feynman diagrams for the process gg HH with the contribution 
from the first top quark KK modes, denoted by tlkkl and t2kkl in the figure. 
Permutations of the external lines are not shown. 



where 

Asm = [{§ - 4m2)Co(s, mf) ~ 2] , (16) 
Akk = 2mt^mt,„ sin [(s - 4m^_„)Co(s, to?_„) - 2] , (17) 

n 

and, as usual 

Co(s,m2)= -| [arcsin(l/V7)]^ for r > 1 

2 

for r < 1 (18) 

with T = 4m^/s. Notice that the triangle contribution is independent of the 
Higgs boson mass. The factor of 2 in eq. (fT7|) is due to the presence of 2 top 
KK excitations for each level. Figure [3] shows the analytical result for difference 
in the triangle contribution only, showing the rapid convergence of the result 
and its agreement with the numerical computation shown in Figured 

The box contribution is more difficult to analyze due to the fact that many 
Passarino-Veltman integrals with different arguments appear in the result. This 
is the reason of the more complicated behavior of the box contribution depicted 
in Figure!^ We notice that there is a strong interference between the triangle 
and box contributions. The final result shows large deviations both enhancing 
and suppressing the cross section, depending on the Higgs boson mass. In the 
case of a Higgs boson lighter than AIh ~ 200 GeV, the partonic gluon fusion 
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Figure 2: Deviations from SM arising separately from triangle and box con- 
tributions, together with the total deviation, as a function of the Higgs boson 
mass. Center-of-mass energy is fixed at %/s = 1 TeV and compactification scale 
l/i? = 500 GeV. 

cross section can be enhanced by more than 40%. These deviations increase 
with partonic center-of-mass energy. 

We present in Figure H] the deviations from the SM result for the partonic 
gluon fusion Higgs pair production for different values of the compactification 
scale for a fixed value of the center-of-mass energy at ^/l = 1 TeV. As expected, 
for larger values of 1/i? the KK modes get heavier and the deviations from the 
SM rapidly decreases. 

The total Higgs pair production cross section at the LHC is computed in the 
standard way by convoluting the partonic cross section with the gluon distribu- 
tion function. We used the Mathematica package implementation of the parton 
distribution functions of [13j with factorization and renormalization scales given 
by — s, as{Mz) = 0.118 and rrit = 175 GeV. In Figures and [S] we com- 
pare the SM result as a function of the Higgs mass with the UED results with 
compactification scales of 1/i? = 500, 700 and 1000 GeV. Differences as large as 
±40% can arise in these models. 

Electroweak precision data puts bounds on the compactification scale as a 
function of the Higgs mass and these constraints decrease with increasing Higgs 
mass [3] . These constraints are sensitive to the top mass allowing — 600 GeV 
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Figure 3: Analytical calculation of deviations from SM arising only from the 
triangle as a function of the KK level. Center-of-mass energy is fixed at Vs — 1 
TeV and compactification scale 1 /R = 500 GeV. 

for niH = 115 GeV and rrit = 173 GeV , which increases in 23% the SM cross 
section. Bounds coming from & — > 57 process [12] implies a compactification 
scale as low as 600 GeV independent of the Higgs mass. The cross section is 
increased by 16% for a light Higgs with mass 120 GeV and = 700 GeV. 

4 Conclusions 

In this paper we studied the effects of UED in the gluon fusion Higgs pair 
production cross section. We implemented the contributions of the top KK 
excitations for the triangle and box diagrams and showed that the partonic 
cross section shows large deviations both enhancing and suppressing the cross 
section, depending on the Higgs boson mass. The total gluon fusion Higgs pair 
production cross section at the LHC can be modified by up to 23% when bounds 
from precision measurements are taken into account. These effects are rapidly 
reduced for larger values of the compactification scale. 
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Figure 4: Deviations from SM gluon fusion Higgs pair production arising from 
top KK modes as a function of the Higgs boson mass for different compactifi- 
cation scales. Center-of-mass energy is fixed at ^/l = 1 TeV. 
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Figure 6: Deviations from SM of Higgs pair production cross section via gluon 
fusion at the LHC as a function of the Higgs mass for values of the compactifi- 
cation scales of 1/i? = 500, 700 and 1000 GeV. 
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